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Pasadena,  California 

Abstract 

It  is  well-known  in  sonar  work  that  the  pulse  form  of  a  direct 
echo  from  a  target  bears  little  relation  to  the  form  of  the  original  signal. 
This  is  true  even  for  regularly  shaped  bodies,  such  as  a  sphere.  In  this 
paper,  the  case  of  a  homogeneous  elastic  sphere  in  water  is  examined 
theoretically  and  it  is  shown  in  comparison  with  experimental  results, 
that  the  observed  effects  originate  from  vibrations  induced  in  the  sphere 
by  the  incident  sound.  Calculated  results  are  presented  for  a  variety  of 
solid  materials  and  it  seems  that  echo  forms  could  possibly  provide  in¬ 
formation  about  the  size  and  constitution  of  a  sonar  target. 
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1 .  Introduction 

It  is  well-known  in  sonar  work  that  the  pulse  form  of  the  direct 
echo  returned  by  a  stationary  insonified  target  in  water  is  usually  quite 
different  from  that  of  the  original  signal  sent  out  by  the  transducer.  This 
effect  can  be  observed  even  when  the  target  has  a  regular  shape  as  in  the 
case  of  a  sphere.  In  the  experiments  which  have  been  made,  the  incident 
sound  has  consisted  of  single  frequency,  constant  amplitude  pulses  of 
various  lengths,  and  the  echo  pulse  generally  appears  in  the  form  of 
multiple  echoes  of  the  original  pulse;  i.e.  compared  to  the  original  pulse, 
the  echo  is  generally  longer  and  subject  to  amplitude  modulation.  Pre¬ 
sumably  there  are  also  differences  in  frequency  content,  but  there  does 
not  appear  to  be  any  quantitative  data  available  on  the  subject. 

If  the  body  has  an  irregular  shape  it  is  possible  to  suppose  that 

this  effect  is  due  to  echoes  returned  by  the  individual  irregularities. 

However  in  the  case  of  regularly  shaped  bodies  with  no  abrupt  changes  in 

curvature,  such  an  explanation  cannot  be  used.  In  this  event  it  would 

seem  reasonable  to  suppose  that  the  distortion  in  the  echo  is  caused 

either  by  diffraction  or  by  vibrations  occurring  within  the  solid  material 

of  the  target  or  by  both.  The  frequencies  used  in  sonar  usually  preclude 

the  influence  of  diffraction,  so  that  the  observed  effects  would  appear  to 

be  due  mainly  to  vibrations  in  the  solid.  Since  the  density  of  any  solid 

does  not  differ  from  that  of  water  by  much  more  than  a  factor  of  eight,  it 

seems  quite  possible  for  the  incident  sound  to  cause  vibrations  in  the 

solid  material  of  the  target.  In  air  the  corresponding  density  ratio  would 

4 

be  of  the  order  of  10  so  that  a  target  would  react  more  like  a  rigid  body, 
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with  a  consequent  diminution  in  echo  distortion. 

It  is  the  purpose  of  this  paper  to  test  the  validity  of  this  hypothesis 
in  the  case  of  a  homogeneous  solid  sphere  supporting  shear  and  com- 
pressional  waves.  Suitable  experimental  data^^has  recently  become 
available  and  this  is  compared  with  calculated  results  based  on  known 

formal  solutions These  results  were  obtained  using  a  high  speed 

[4]  [5] 

computer.  Previous  calculations  have  been  made  for  fluid1-  1  and  rigid1-  J 
spheres . 

2.  Formulation  of  the  Problem 

The  coordinate  system  for  the  sphere  is  shown  in  Fig.  1  where 
the  relationship  between  the  cartesian  and  spherical  polar  coordinates  is 

x  =  r  sin  9  cos  <f>  , 
y  =  r  sin  9  sin  <f)  , 

z  =  r  cos  9  .  (1) 

The  sphere  is  assumed  to  consist  of  solid  isotropic  material  supporting 
both  compressional  and  shear  waves.  The  displacement  vector  u  can 
be  expressed  using  the  vector  and  scalar  potentials  A  and  ij*  as 
follows  ^  ^ 


[lj  L. D.  Hampton  and  C. M.  McKinney;  J .  A. S.  A.  ,  33,  5,1961,  p.  664. 

[2]  J.  J.  Faran;  J .  A. S.  A. ,  23,  4,  1951,  p.  405. 

[3]  P.  M.  Morse  and'H.  Feshbach;  "Methods  of  Theoretical  Physics" 

(McGraw-Hill  Book  Company,  New  York,  1953)  Vol.  II,  p.  1483. 

[4]  V.  C.  Anderson;  J.A.S.A.  22,  4,  1950,  p.  426. 

[5]  H.  Stenzel;  Leitfaden  zur  Berechnung  von  Schallvorgagen,  Julius 

Springer,  Berlin,  1939. 

[6]  Morse  and  Feshbach  p.  142 
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where 


u  =  -Vi{;  +  V*  A  , 


V*i}»=  (1  /c^2)  02  i|i/9t2  , 

V2A  =  (1/c  2)  82  A/8t2  , 
2 


(2) 

(3) 

(4) 


describe  the  motion  of  the  compressional  and  shear  waves  respectively, 


c^  and  c^  are  the  compressional  and  shear  wave  velocities  defined  by 


=  [E(l-cr)/ (l+o)(l~2cr)]2 

c  =  [E/2p  (l+or)p  ,  (5) 

2  1 


where  E,  ,  o  are  the  Young's  modulus,  density  and  Poisson's  ratio 
of  the  solid  material  of  the  sphere. 

Outside  the  sphere  there  is  a  limitless  fluid  of  density  p  and 
sound  velocity  c  in  which  there  is  a  continuous  train  of  waves  emanating 
from  a  point  source  situated  on  the  z  axis  at  r  =  rQ>  0  =  it.  The  time 
dependence  of  these  waves  is  of  the  form  exp(-iwt)  from  which  the  wave 
number  k  in  the  fluid  is  obtained  by  means  of  the  relation 


k  =  w/  c  =  2ir/X 


where  X  is  the  wave  length.  Similar  relations 


k  =  «/  c  ;  k  =  <o/  c 

112  2 


hold  for  the  compressional  and  shear  waves  in  the  solid.  The  waves 

[  7 1 

emanating  from  the  point  source  can  be  expressed1-  J 


[7]  ibid  p.  1466 
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p^  =  Pq  exp(ikD)/D 
oo 

=  ikPo  ^  (2n+l)(-l)nPn(cos  0)jn(kr)hn(krQ) 
n=o 

where 


0  <  r  <  r 

o 


(6) 


i 

D  =  [r*  +  2rr  cos  0  +  r2l2 
o  o 


and  the  Pn  are  Legendre  polynomials,  and  the  j  ,  hn  are  spherical 
r  3 1 

Bessel  functions1  .  Plane  waves  incident  on  the  sphere  are  obtained  by 
making  rQ  go  to  infinity.  Using  the  two  limits 

exp  (ikD)/D — ►exp  (ik  rQ)  exp  (ik  r  cos  0)/rQ 

hn(krQ)  i“(n+1)  exp  (ik  rQ)/kro  (7) 


and  removing  the  common  factor  exp(ik  rQ)/ rQ  gives 

p.  =  PQ  exp  (ikr  cos  0) 
oo 

=  PQ  ^  (2n+a)inPn(cos0)jn(kr)  .  (8) 

n=o 

The  above  waves,  p.,  incident  on  the  sphere  result  in  scattered 
waves  in  the  fluid  which  will  be  of  the  form 

oo 

Ps  =  Po  I  cnhn(kr)Pn(cos  0)  (9) 

n=o 


[8]  ibid  p.  1325  and  p.  1573 
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where  the  coefficients  c^  have  to  be  determined  from  the  boundary  con¬ 
ditions  at  the  surface  of  the  sphere.  The  appropriate  conditions  are: 

(a)  the  pressure  in  the  fluid  is  equal  to  the  normal 
component  of  stress  in  the  solid. 

(b)  the  normal  component  of  displacement  of  the 
fluid  is  equal  to  the  normal  component  of  displace¬ 
ment  of  the  solid. 

(c)  the  tangential  components  of  shearing  stress  in 
the  solid  vanish  at  the  surface. 

These  conditions  can  be  expressed^ 


p.  +  p  +  2p  c 
*s  2 


f9l 

J  in  the  relations 

(  8u  ' 

a 

divu  +  “UF 

(1-2*) 

=  0 


u.  +  u  =  u 
i,  r  s,  r  r 


and 


3u0 

8r 


u 


e 


,  1  r 
+  T  W 


=  o 


(10) 

(ii) 


(12) 


evaluated  on  the  surface  of  the  sphere  at  r  =  a.  The  displacement  of  the 
fluid  in  Eq.  (11)  is  obtained  from  (6)  and  (9)  by  using 


"i,r+us,r  8<P1+Ps,/8r  • 


(13) 


Because  of  the  symmetry  of  the  incident  waves,  the  component  of  dis¬ 
placement  u^  is  taken  to  be  zero  and  the  <f>  component  of  the  shearing 

stress  can  also  be  neglected.  The  only  non-zero  component  of  the  vector 


[9]  ibidp.  142  and  p.  116. 
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potential  will  be  so  that  the  potentials  become 

oo 

*=  I  Vn<V)Pn(c0s9) 

n=o 

oo 

'  I  WV’  Pn<cos  *>  •  (14) 

n=o 

Using  Eq.  (2)  these  expansions  can  be  inserted  into  the  boundary  condi¬ 
tions  (10)  -  (12)  and  the  coefficients  of  each  normal  mode  equated.  The 
coefficients  cn  of  the  scattered  waves  in  the  fluid  given  by  Eq.  (9)  can 
then  be  expressed  as  follows, 

cn  =  k(-l)n(2n+i)hn(kro)  sinr^  expt-ir^)  (15) 

where  the  angle  is  given  by 

tanr^  =  -  [jn(x)Fn-j'n(x)]/[nn(x)Fn-n'n(x)]  (16) 

with 

2(n2+n)jn(xz) 

(n2+n-  2)jn(xz  )+xz2  jn"(x2 ) 
Z(n2+n)[jn(x2)-x2jn'(x2)] 

(n2  +n~  2)  jn(x2  )+X22  jn"(x2 ) 

(17) 

x  =  ka  ;  x  =  k  a  ;  x 

1  1 


F  = 
n 


x 

-4 

p,  ~ 


X  j  '(x  ) 
1  n  i ' 


l2\lT2o^Xl)^n{xl)] 


and 


=  k  a  . 
2  2 
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The  primes  denote  the  derivative  with  respect  to  the  argument.  This 

r  21 

result  was  first  derived  by  FaranL  .  However  there  was  an  error  in 
his  presentation  in  which  the  factor  o/(l-2a)  was  misplaced.  Finally  it 
should  be  noted  that  the  expression  on  the  right  hand  side  of  Eq.(15)  is  of 
the  form  f(x)/(f(x)  +  i  g(x) )  where  f,  g  are  regular  on  the  real  axis. 
Hence  there  are  no  singularities  when  the  argument  is  real,  and  the 
function  can  be  integrated  numerically  in  a  straightforward  manner.  It 
also  follows  that  the  solution  as  presented  is  complete  for  all  frequencies, 
i.e.  ,  the  boundary  conditions  are  fully  satisfied  by  the  shear  and  com- 
pressional  waves  postulated  in  Eq.  (14). 

Certain  limiting  cases  are  of  interest.  If  Fn-”*-0,  the  solution 
would  then  apply  to  scattering  by  a  rigid  immovable  sphere^ This 
would  be  the  case  for  instance  when  the  density  of  the  solid  was  very 
much  greater  than  that  of  the  fluid.  If  F  oo  the  solution  for  scatter¬ 
ing  by  a  free  surface  sphere  is  obtained^  This  corresponds  to  the 
condition  where  the  normal  stress  at  the  surface  of  the  sphere  vanishes, 
which  would  result  for  example  when  the  density  of  the  material  inside 
the  sphere  was  very  much  less  than  that  of  the  fluid. 

From  the  above  it  follows  that  the  echo  returned  by  the  solid 
sphere  to  the  source  is  given  by 


[lOj  P.  M.  Morse;  Vibration  and  Sound  (McGraw-Hill  Book  Co.  Inc. 
New  York  1948)  p.  354 

[ll]  Morse  and  Feshbach  p.  1483 
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P 

Pe  =  ~Z¥~ 

o 


00 

2xq  ^  (2n+l)  sin exp(- i t^)hn*(x0) 


n=o 


exp(-ioot) 


=  2r“  f(x,xQ,x  ,x  )  exp  (-ixT) 

o  i  a 


(18) 


where  xq  =  krQ  =  xR,  and  T  =  ct/a.  When  the  source  is  a  large  distance 
from  the  sphere 


Pe  = 


jV 

2r  2 
o 


oo 

\  ^  (-l)n(2n+l)  sin  exp  (-ir^) 


L  n=o 


exp  [ik(2rQ-ct)  ] 


P  a 
o 

2r  2 
o 


yx.x  .x^)  exp  [ix(2R-r)]  . 


(19) 


Removal  of  a  factor  exp(ikro)/rQ  gives  the  solution  for  incident  plane 
waves.  Equations  (18)  and  (19)  can  then  be  used  to  construct  the  echo  due 
to  a  pressure  pulse  emanating  from  the  source.  Suppose  the  source 
emits  a  pulse  of  form  P-(t).  This  can  be  expressed  in  terms  of  Fourier 
components  as 


P  c  n  oo 

P4(t)  =  - — -  J  g(k)  exp  [ik(D-ct)  ]  dk 


(20) 


( 2ir) ZD  "°° 


where  D  is  as  defined  previously.  The  frequency  spectrum  g(k)  is 
found  by  taking  the  Fourier  transform  of  the  given  pulse,  i.e. , 

— f  P^t)  exp  [-ik(D-ct)]  dt  . 

(2tt)2  ■°° 


g(k)  = 


(21) 
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With  the  new  variable  x  =  ka,  the  reflected  pulse  will  be 

r 


P  c 

p  (T)  =  - - -  \  g(x)  f  (x,x  ,x  ,x  )  exp  (-ixr)  dx  , 

C  i-  v  U  1  i 

2  ,  -  CO 


(22) 


(2-rr)zar 


o 


and  when  the  point  source  moves  to  a  large  distance  from  the  sphere 


P  .  oo 

PJ-7-)  =  - 2^ -  \  g(x)f  (x,x  ,x  )  exp  [ix(2R-T)]  dx 


2(  2tt) 2  r 


2^  2  ''-oo 


(23) 


In  general  the  echo  given  by  (23)  will  differ  in  form  from  that  of 

the  incident  pulse  (20).  Only  for  high  frequencies  in  the  special  cases  of 

f  i  2] 

a  rigid  and  a  free  surface  sphere  will  it  be  the  same.  It  can  be  shown1  J 
that  in  the  former  case  f ^ — ►exp(-2ix)  as  x  becomes  large  while  for  the 
free  surface  sphere  f^^— *■ -exp  (- 2ix) .  Hence  if  the  frequencies  contained 
in  the  pulse  are  in  the  high  frequency  range,  Eq.  (23)  becomes 

P  n  oo 

Pr(T)  =  - - -  \  g(x)  exp  [  ix(2R-  2-  r)]  dx 

2(2Tr)2roZ  "°° 

which  means  that  the  reflected  pulse  has  the  same  form  as  the  emitted 
pulse, but  is  returned  time  2(rQ-a)/c  later.  This  travel  time  indicates 
that  the  sound  is  reflected  from  a  point  source  reflector  at  the  point  on 
the  surface  of  the  sphere  nearest  to  the  source  of  incident  sound.  A 
similar  result  holds  for  the  free  surface  sphere  except  that  the  pulse  is 
inverted. 


[12]  ibid  p.  1554. 
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3.  The  Steady  State  Solutions 

The  steady  state  solutions  given  by  the  functions  f  in  Eqs.  (18), 
(19)  were  determined  for  a  certain  number  of  cases,  the  calculations  being 
performed  on  a  high  speed  computer.  The  results  are  shown  Figs.  2-13. 

The  first  results  obtained  were  for  the  special  cases  of  the 
rigid^*^  and  free  surface  ^  ^spheres  for  a  distant  point  source  of  con¬ 
tinuous  waves.  These  are  shown  in  Figs.  2(a),  (b).  The  argument  of  the 

function  f  is  presented  divided  by  the  variable  x  =  ka.  The  results  for 

r  5 1 

the  rigid  sphere  are  in  agreement  with  those  of  StenzelL  .  For  low 
frequencies  the  function  is  given  by  the  initial  terms  in  the  series  expan¬ 
sion  which  in  the  limit  as  x  tends  to  zero  are 


and 


2i 

x 


THIS) 


3jj'(x) 

TyTST 


2i  jo<x) 

X  hQ(x) 


2(1 -ix)  , 


for  the  rigid  and  free  surface  spheres  respectively.  For  the  rigid  sphere 
this  represents  the  well-known  condition  of  Rayleigh  scattering  where  the 
scattered  intensity  is  proportional  to  the  fourth  power  of  the  frequency. 
For  the  free  surface  sphere  the  results  are  quite  different.  Not  only 
does  the  scattered  intensity  reach  maximum  values  at  low  frequencies, 
but  the  scattering  is  uniform  in  all  directions.  For  high  frequencies  both 
solutions  tend  to  form  exp(-2ix),  the  free  surface  solution  converging 
more  rapidly  than  that  of  the  rigid  sphere.  In  the  previous  section,  it 
was  shown  that  this  indicates  that  at  high  frequencies  the  sound  is 
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mainly  from  a  small  area  on  the  surface  opposite  the  source  and  this 
would  agree  with  physical  intuition.  At  low  frequencies  the  echo  appears 
to  come  from  the  center  of  the  rigid  sphere  and  from  a  half  radius  posi¬ 
tion  in  the  free  surface  spliere.  As  the  frequency  increases  the  apparent 
origin  of  the  echo  moves  gradually  towards  the  region  on  the  surface 
opposite  the  source.  This  is  shown  in  Fig.  2(b),  where  the  phase  of  f 
divided  by  x  =  ka  represents  distance  along  a  radius  inside  the  sphere. 

In  the  case  of  the  rigid  sphere  these  results  can  be  readily  understood  by 
supposing  that  low  frequency  waves  are  intercepted  by  the  entire  cross 
section  whereas  high  frequency  waves  behave  as  in  geometrical  optics  and 
form  a  "bright  spot"  reflector  on  the  surface  opposite  the  source.  With 
the  free  surface  sphere  the  results  at  low  frequencies  are  not  so  readily 
explainable  except  that  as  expected  they  differ  from  those  for  the  rigid 
sphere. 

The  main  body  of  results  were  derived  for  solid  spheres  support¬ 
ing  internal  shear  and  compressional  waves.  The  properties  of  the 
materials  considered  are  given  in  Table  I.  The  fluid  outside  the  sphere 
was  assumed  to  be  water  of  density  1  gm/cc  and  compressional  velocity 
1 , 41 0  m/sec. 

[  2 

As  an  initial  test  of  the  programs,  the  results  obtained  by  FaranL 
were  recalculated.  Since  these  were  for  a  -  y  no  error  could  result 
from  the  misplacing  of  the  factor  o7(1-2ct)  mentioned  in  the  previous 
section  since  this  factor  is  unity.  Good  agreement  was  found. 

Some  of  the  results  for  the  materials  tested  in  Table  I  are  given 
in  detail  in  Figs.  3-11.  As  before  these  are  for  a  distant  source. 
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TABLE  I^-13^ 


Material 

Density 

(gm/cc) 

Poisson's 
Ratio  a 

Co  mp  r  e  s  s  i  onal 
Velocity 

C^  (m/sec) 

Shear  Velocity 

C  (m/sec) 

2 

Beryllium 

1.87 

0.  05 

12,890 

8,880 

Fused  Silica 

2.  20 

0.17 

5,968 

3,764 

Heavy  Silicate, 
Flint  Glass 

3.88 

0.  224 

3,980 

2,  380 

Armco  Iron 

7.  70 

0.  29 

5,960 

3,  240 

Monel  Metal 

8.90 

0.  327 

5,  350 

2,720 

Aluminum 

2.70 

0.  355 

6,420 

3,040 

Yellow  Brass 

8.  60 

0.  374 

4,700 

2,110 

Lucite 

1.18 

0.40 

2,  680 

1,100 

Lead 

11. 34 

0.43 

1,960 

690 

IceM 

0.917 

0.  336 

2,743 

1,433 

Generally  the  range  of  frequency  was  for  values  of  ka  up  to  30,  but  for 
Armco  Iron,  ice  and  lucite  the  range  extended  to  ka  =  60,  20,  10  respec¬ 
tively.  In  addition  to  the  pressure  amplitude  the  phase  variation  is  given 
for  Armco  Iron,  and  aluminum.  In  all  cases  the  results  begin  at  low 
frequencies  as  though  the  solid  were  a  rigid  body,  changing  in  general 

[13]  American  Inst,  of  Physics  Handbook  (McGraw-Hill  Book  Company, 

Inc.  New  York,  1957). 

[14]  D.  L.  Anderson;  Trans,  of  Engin.  Inst,  of  Canada,  Vol.  2,  pp.  1 1 6-1  22, 

1958. 
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into  a  fairly  regular  series  of  peaks  and  minima  as  the  frequency  in¬ 
creases.  With  a  rigid,  incompressible  material  such  as  beryllium,  the 
change  from  the  rigid  body  solution  is  not  very  great.  However,  as  the 
material  becomes  more  compressible  and  pliant,  the  resonances  tend  to 
become  more  pronounced  and  more  closely  spaced.  In  the  case  of  lucite 
and  ice  the  resonances  have  become  quite  sharp  and  close  together.  This 
general  trend  was  investigated  by  considering  an  average  frequency  in¬ 
terval  between  minima  or  between  resonance  peaks  for  each  material. 

The  results  are  shown  in  Fig.  12  plotted  against  the  shear  velocity  c^. 

Parameters  other  than  c  were  also  considered  such  as  Poisson's  ratio 

2 

c r,  but  with  these  the  scatter  of  points  was  much  greater.  It  appears 

therefore  that  this  feature  is  most  strongly  dependent  on  the  behavior 

of  shear  waves  in  the  material.  The  successive  peaks  and  minima  which 

occur  in  the  direct  echo  for  a  continuous  frequency  were  shown  by 

[2  1 

Faran  J  to  be  due  to  strong  lobes  of  backscattered  radiation  forming  and 
then  splitting  again  into  side  lobes  scattering  in  other  directions.  It 
seems  therefore  that  shear  waves  play  an  important  role  in  this  process. 
Other  general  trends  can  be  observed  in  Figs. 3  -  11.  Peaks  split  off  and 
merge  as  the  properties  of  the  material  change,  and  as  before  it  appears 
that  the  shear  wave  velocity  is  the  dominant  parameter.  The  density  of 
the  material  seems  to  have  a  relatively  unimportant  effect  in  water.  As 
the  frequency  increases  there  appears  to  be  no  tendency  towards  some 
constant  limit  as  in  the  cases  of  the  free  surface  and  rigid  spheres.  The 
peaks  seem  to  recur,  but  in  an  increasingly  ragged  form. 

In  the  phase  variations  shown  in  Figs.  6,  8  jumps  in  phase 
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occur  at  frequencies  corresponding  to  minima  in  the  pressure  amplitude. 
Although  the  parameter  (-arg  f ^ka)  varies  continuously  from  zero  as  in 
the  case  of  the  rigid  sphere,  it  does  not  appear  that  any  direct  physical 
interpretation  should  be  applied  to  it  concerning  the  apparent  source  of 
the  echo.  The  representation  should  be  regarded  only  as  a  convenient 
way  of  presenting  the  phase  as  a  continuously  varying  function. 

In  order  to  determine  the  effect  of  distance  of  the  sound  source 
from  the  center  of  the  sphere  the  function  f  in  Eq.  (18)  was  evaluated 
for  a  point  source  at  various  distances  from  a  rigid,  a  free  surface,  and 
a  brass  sphere.  The  pressure  amplitude  j  f  |  multiplied  by  R  is  shown 
in  Figs.  1  3  -  15,  allowing  a  ready  comparison  with  the  solutions  for  a 
distant  source  |  f^J.  In  general  it  appears  that  f  represents  a  satis¬ 
factory  solution  when  the  sound  source  is  situated  more  than  10  radii 
from  the  center  of  the  sphere. 

4.  Echo  Pulse  Forms 

The  most  obvious  general  feature  in  the  steady  state  solution  for 
ordinary  metals  is  the  succession  of  peaks  and  minima  in  the  pressure 
amplitude  and  it  is  of  interest  to  determine  how  this  affects  the  pulse 
form  of  the  echo  when  the  steady  state  solutions  are  used  in  the  integral 
expression  (23)  for  a  distance  source.  The  incident  pulse  form  could  be 
chosen  arbitrarily.  However  in  practice  the  incident  sound  is  generally 
produced  by  making  a  transducer  resonate  over  several  cycles  at  a 
particular  frequency.  Mathematically  the  pressure  variation  which 
results  at  a  point  in  the  fluid  can  be  represented  as  follows, 
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P.(t)  =  0 

=  exp  (-i  WQt) 

=  0 

where  coq  is  the  angular  frequency  of  the  transducer  at  resonance,  and 
2  A  t  is  the  duration  of  the  pulse.  The  frequency  spectrum  g(w)  is 
given  by  the  transform 

1  f At 

g(ai)  =  -y-  \  exp  [i(w-aj  )t]  dt 

^  J-At 

2 

=  —  sin  [(w-co  )  At]  /  (w-w  ) 

TT  O  J  o' 

which  in  the  non-dimensionalized  system  of  Eq.  (23)  becomes 

g(x)  =  ~  sin  [  (x  -  xq  )  At]  /  (x  -  xq)  (25) 

where  xq  =  fc^a/c  and  is  referred  to  as  the  dominant  frequency.  By  use 
of  Eq.  (25)  and  the  previously  derived  values  of  the  function  f^  it  is  then 
possible  to  obtain  the  pulse  form  of  the  echo  by  numerically  integrating 
Eq.  (23). 

The  nature  of  the  function  g  in  Eq.  (25)  is  shown  in  Figs.  16  - 
1  8  for  different  pulse  lengths  A  T.  The  height  of  the  main  peak  occurring 
at  the  dominant  frequency  is  equal  to  At  and  its  "spread"  varies  inverse¬ 
ly  with  A  T.  If  the  function  |  f  |  is  momentarily  idealized  as  consisting 
of  a  series  of  similar,  equally  spaced  peaks,  it  would  appear  that  the 
form  of  the  echo  depends  mainly  on  the  pulse  length  and  on  the  location 
of  the  dominant  frequency  relative  to  the  maxima  and  minima  of  the  |f  |. 


t  <  -  At 

-  At  <  t  <  At  (24) 

t  >At 
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Two  extreme  cases  would  then  arise  depending  on  whether  the  dominant 
frequency  coincided  with  a  maximum  or  with  a  minimum  of  |f  |. 

Using  the  data  for  Armco  iron  as  shown  in  Figs.  5,  6  several 
echo  pulse  forms  were  computed  for  different  lengths  of  the  incident  pulse 
and  for  dominant  frequencies  corresponding  to  values  of  x  or  ka  at 
24.5  and  25.5.  The  former  frequency  occurs  at  a  peak  of  the  pressure 
amplitude  and  the  latter  at  a  minimum.  The  range  of  integration  over  ka 
for  the  longer  pulses  was  from  15  to  35,  while  for  the  short  5  cycle  pulse, 
it  extended  from  10  to  40.  The  incident  pulse  did  not  therefore  have  a 
perfectly  rectangular  form.  However  in  comparison  to  experimental 
pulse  forms,  it  could  be  considered  a  satisfactory  approximation.  In 
addition  the  irregularities  introduced  by  restricting  the  range  of  integra¬ 
tion  facilitated  the  recognition  of  certain  features  of  the  incident  pulse  in 
the  echo.  The  calculated  echoes  are  shown  in  Figs.  16  -  18.  The  time 
scale  for  the  incident  pulse  is  chosen  with  respect  to  the  time  of  arrival 
of  the  mid-point  of  the  pulse  at  the  center  of  the  sphere,  whereas  the 
scale  for  the  echo  is  chosen  with  respect  to  the  time  of  arrival  back  at 
the  source.  All  pulses  are  shown  travelling  from  right  to  left.  It  can  be 
seen  from  these  figures  that  the  leading  edge  of  the  echo  precedes  that  of 
the  incident  pulse  by  a  time  difference  of  2  in  each  case.  In  addition 
the  leading  edge  of  the  echo  is  of  the  same  form  as  the  leading  edge  of 
the  incident  pulse.  These  features  indicate  that  the  first  part  of  the  echo 
consists  of  a  rigid  body  reflection  from  the  region  of  the  surface  of  the 
sphere  adjacent  to  the  sound  source.  The  subsequent  parts  of  the  echo 
are  affected  by  the  vibrations  of  the  sphere.  In  the  case  of  the  five  cycle 
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pulse,  the  first  echo  is  of  identical  form  to  the  incident  pulse,  while  the 
second  echo  is  also  of  the  same  form,  but  inverted.  Subsequent  echoes 
diminish  in  amplitude  and  lose  the  characteristic  features  of  the  incident 
pulse.  Whether  the  dominant  frequency  occurs  at  a  minimum  or  a  maxi¬ 
mum  of  the  function  |  f  |  does  not  appear  to  make  much  difference  to 
the  form  of  the  echo  for  the  short  five  cycle  pulse,  but  obviously  it  is 
important  when  the  pulse  is  longer.  The  reason  can  be  seen  from  the 
frequency  spectra  shown  in  Figs.  16-18.  A  change  in  ka  of  the  order 
of  1  in  the  spectrum  for  the  short  pulse  will  not  greatly  affect  the 
integral  (23);  however,  this  is  not  the  case  for  the  longer  pulses.'  The 
differences  in  form  of  the  echoes  shown  in  Figs.  17,  18  are  in  fact  quite 
distinctive.  It  seems  moreover  that  the  changes  which  occur  when  the 
dominant  frequency  is  moved  from  a  maximum  of  |  f ^  |  to  a  minimum, 
are  characteristic  of  any  ordinary  metal.  Figure  19  shows  the  echoes 
from  a  sphere  of  aluminum  for  the  same  type  of  incident  pulse,  where 
these  results  were  obtained  experimentally.  With  allowance  for  a  change 
of  scale  and  such  effects  as  the  response  of  the  transducer,  the  differences 
in  the  echo  resulting  from  a  change  in  the  dominant  frequency  are  closely 
related  to  those  shown  in  Figs.  16  -  18.  The  sphere  in  this  case  had  a 
diameter  of  5  inches  and  the  change  in  dominant  frequency  was  from 
120  kc/sec  to  123.5  kc/sec.  This  is  equivalent  to  a  change  in  ka  of 
about  1,  which  according  to  Fig.  12  is  the  approximate  distance  be¬ 
tween  a  peak  and  a  minimum  of  the  function  jf  |  for  most  metals  in¬ 
cluding  aluminum.  Using  the  constants  given  above  for  water  and 
aluminum,  it  is  found  that  120  kc/sec  does  not  in  fact  coincide  with  a 
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120  kc/sec  to  123.5  kc/sec.  This  is  equivalent  to  a  change  in  ka  of 
about  1,  which  according  to  Fig.  12  is  the  approximate  distance  be¬ 
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aluminum,  it  is  found  that  120  kc/sec  does  not  in  fact  coincide  with  a 
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peak  of  the  steady  state  reflection  function  |f  j.  However  this  is  not 
surprising  since  the  values  used  referred  to  rolled  aluminum.  In  addi¬ 
tion  if  the  frequency  is  to  be  expressed  in  terms  of  ka  with  any  accuracy, 
it  would  be  necessary  to  know  the  velocity  of  sound  in  water  under  the 
conditions  of  the  experiment,  and  also  the  diameter  of  the  sphere,  to 
within  less  than  1%.  Echoes  were  calculated  for  rolled  aluminum  at 
values  of  ka  equal  to  34.  6  and  35.  6  corresponding  to  frequencies  of 
122.  3  kc/sec  and  125.8  kc/sec,  these  values  occurring  at  a  maximum 
and  a  minimum  respectively  of  the  reflection  function  |  f  |  .  The  pulse 
lengths  A  T  were  the  same  as  those  in  Figs.  16-18.  Similarly  echoes 
were  calculated  for  a  brass  sphere  for  frequencies  at  ka  =  20.  2  and 
21.0.  These  echoes  were  found  to  have  the  same  features.  The  leading 
edge  was  a  rigid  body  reflection  and  the  same  kind  of  transition  in  pulse 
form  occurred  when  the  frequency  and  the  pulse  length  were  varied. 

With  yellow  brass  the  secondary  echoes  in  the  multiple  echo  forms  had  a 
bigger  amplitude  than  the  primary  echo. 

5.  Discussion 

Although  this  paper  represents  only  a  preliminary  study,  it  may 
be  worthwhile  to  consider  the  significance  of  the  results  in  relation  to 
the  problem  of  using  sonar  echoes  to  obtain  information  about  a  target. 

In  the  first  place,  it  would  seem  that  solid  materials  could  be 
divided  roughly  into  two  groups,  metallic  flint-like  substances  and  sub¬ 
stances  which  are  fairly  pliant.  This  can  be  seen  from  the  steady  state 
solutions  where  the  former  is  characterized  by  a  succession  of  peaks 
and  minima  roughly  the  same  distance  apart,  while  the  latter  has  sharper, 
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stronger  peaks  more  closely  spaced.  Although  all  the  echo  forms  which 
were  calculated  belong  only  to  the  first  type,  it  is  evident  from  the  steady 
state  solutions  that  there  would  be  a  difference  in  the  general  nature  of 
echoes  between  the  two  groups.  Hence  there  would  exist  the  possibility 
of  distinguishing  for  instance  between  a  bare  rock  and  a  large  fish. 

Secondly,  if  the  sonar  target  is  known  to  be  a  homogeneous 
metallic  sphere,  then  it  is  possible  to  determine  its  approximate  radius 
by  using  data  of  the  type  shown  in  Fig.  19.  The  features  of  the  transition 
between  a  peak  and  a  minimum  of  the  steady  state  reflection  function 
|  f  |  for  the  long  incident  pulses  is  characteristic  of  most  ordinary 
metals,  as  shown  in  the  previous  section.  The  transition  is  accomplished 
during  a  change  in  ka  of  the  order  of  1.  Hence  given  an  actual  change 
in  frequency  in  cycles  per  second,  it  is  then  possible  to  determine  the 
radius  a  of  the  sphere.  For  example  the  transition  in  pulse  form 

shown  in  Fig.  1  9  is  achieved  through  a  frequency  change  of  3.  5  kc/sec 
corresponding  to  a  change  in  ka  of  about  1  which  therefore  makes  the 
radius  of  the  sphere  approximately  2.5  inches.  In  general  however  such 
an  estimate  would  not  be  quite  so  accurate.  It  might  also  be  possible  to 
estimate  the  size  of  the  sphere  by  varying  the  pulse  length  rather  than 
the  frequency,  particularly  if  the  dominant  frequency  occurs  at  a  mini¬ 
mum  of  the  reflection  function  If  I.  However  it  would  presumably  be 

'  oo  ’ 

necessary  to  use  both. 

Finally,  it  has  been  shown^^4hat  there  are  significant  differ¬ 
ences  in  the  steady  state  reflection  function  f  for  rigid  bodies  of 


[15]  R.  Hickling;  J.A.S.A.  ,  30,  2,  137-139,  Feb.  1958. 
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different  shapes.  Although  these  effects  would  be  rendered  more  com¬ 
plicated  by  allowing  for  the  vibrations  of  the  solid  material,  it  may  be 
possible  to  use  them  to  derive  some  information  about  the  shape  of  a 
sonar  target. 


Fig.  1  Coordinate  System 


Fig.  2(a)  The  pressure  amplitude  as  a  function  of  frequency  of  the 
echo  returned  by  rigid  and  free  surface  spheres  to  a 
distance  source  of  continuous  waves. 


Fig.  2(b)  The  phase  of  the  echo  from  rigid  and  free  surface 
spheres  as  a  function  of  frequency. 


Fig.  3  The  pressure  amplitude  as  a  function  of  frequency  of  the 
echo  returned  by  a  beryllium  sphere  to  a  distance  source 
of  continuous  waves. 
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Fig.  4 


The  pressure  amplitude  as  a  function  of  frequency  of  the 
echo  returned  by  a  fused  silica  sphere  to  a  distant  source 
of  continuous  waves. 
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rig.  5  The  pressure  amplitude  as  a  function  of  frequency  of 
the  echo  returned  by  a  sphere  of  Armco  iron  to  a 
distant  source  of  continuous  waves. 


Fig.  6  The  phase  of  the  echo  from  the  Armco  iron  sphere  as 
a  function  of  frequency. 
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The  phase  of  the  echo  from  the  Armco  iron  sphere  as 
a  function  of  frequency. 


Fig.  7  The  pressure  amplitude  as  a  function  of  frequency  of 
the  echo  returned  by  an  aluminum  sphere  to  a  distant 
source  of  continuous  waves. 


Fig.  8 


The  phase  of  the  echo  from  the  aluminum  sphere  as 
a  function  of  frequency. 
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Fig.  9  The  pressure  amplitude  as  a  function  of  frequency  of  the 
echo  returnedby  a  sphere  of  yellow  brass  to  a  distant 
source  of  continuous  waves. 


FREQUENCY  INTERVAL,  ka 


SHEAR  VELOCITY,  C2  ,  M/SEC 


Average  interval  in  frequency  between  minima  or 
between  peaks  in  the  pressure  amplitude  of  the  echo 
for  different  materials  as  a  function  of  the  shear 
velocity  of  the  material. 


Fig.  12 


0 


5 

FREQUENCY,  ka 


10 


Fig.  13  The  pressure  amplitude  as  a  function  of  frequency  of 
the  echo  returned  by  a  rigid  sphere  to  a  point  source 
of  continuous  waves  distance  Ra  from  the  center  of 
the  sphere. 
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Fig.  14  The  pressure  amplitude  as  a  function  of  frequency  of 
the  echo  returned  by  a  freee  surface  sphere  to  a  point 
source  of  continuous  waves  distance  Ra  from  the 
center  of  the  sphere. 
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Fig.  i 


The  pressure  amplitude  as  a  function  of  frequency  of 
the  echo  returned  by  a  brass  sphere  to  a  point  source 
of  continuous  waves  distance  Ra  from  the  center 
of  the  sphere. 


FREQUENCY  SPECTRUM  -  DOMINANT  FREQUENCY  ka-24.5 


Fig.  16  Pulse  forms  of  echoes  returned  by  an  Armco  iron 
sphere  for  dominant  frequencies  at  ka  =  24.5  and 
<s5.  5  with  an  incident  pulse  of  five  cycles. 


FREQUENCY  SPECTRUM -DOMINANT  FREQUENCY  ka  -  24.5 


Fig.  19  Experimentally  determined  echo  forms  for  a  5" 
diameter  aluminum  sphere  for  various  incident 
pulse  lengths.  The  sweep  is  150  |U  sec/ cm. 
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